Role of uniform horizontal magnetic field on convective flow 
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The effect of uniform magnetic field applied along a fixed horizontal direction in Rayleigh-Benard 
convection in low-Prandtl-number fluids has been studied using a low dimensional model. The 
model shows the onset of convection (primary instability) in the form of two dimensional stationary 
rolls in the absence of magnetic field, when the Rayleigh number R is raised above a critical value 
R c . The flow becomes three dimensional at slightly higher values of Rayleigh number via wavy 
instability. These wavy rolls become chaotic for slightly higher values of R in low-Prandtl-number 
(P r ) fluids. A uniform magnetic field along horizontal plane strongly affects all kinds of convective 
flows observed at higher values of R in its absence. As the magnetic field is raised above certain 
value, it orients the convective rolls in its own direction. Although the horizontal magnetic field 
does not change the threshold for the primary instability, it affects the threshold for secondary 
(wavy) instability. It inhibits the onset of wavy instability. The critical Rayleigh number R (Q,P r ) 
at the onset of wavy instability, which depends on Chandrasekhar's number Q and P r , increases 
monotonically with Q for a fixed value of P r . The dimensionless number R (Q, P r )/(R c QPr) scales 
with Q as Q" 1 . A stronger magnetic field suppresses chaos and makes the flow two dimensional 
with roll pattern aligned along its direction. 

PACS numbers: 47.20.Bp, 47.20.-k, 47.52.+j, 47.35.Tv, 47.65.-d 



I. INTRODUCTION 

Convective flows in low-Prandtl-number fluids in the 
presence of magnetic field have been studied for many 
years due to its importance in geophysical and astrophys- 
ical problems p], [2j] . Generally this flow is called magneto 
convection. There are also industrial applications of this 
kind of flow in crystal growth 0] and in fusion reactor 
as heat exchanger [1] . There has been extensive theoreti- 
cal and numerical studies of thermal convection in fluids 
in presence of external magnetic field p], [5l4Tl| . These 
studies reveal the stabilizing effect of the magnetic field 
on convective flow. In several experiments [12l-ll4j it have 
been shown that magnetic field affects the convective flow 
strongly. 

Fauve et. al. [H, [l6[ studied the effect of horizontal 
magnetic field, both in the longitudinal and the trans- 
verse directions, on wavy roll instability in a Rayleigh- 
Benard experiment with mercury. They found that the 
horizontal magnetic field inhibited the oscillatory insta- 
bility and made the convection two dimensional. Libch- 
aber et. al. pjj and Hof et. al. pl| considered the ef- 
fect of external horizontal magnetic field on low-Prandtl- 
number thermal convection and found that the magnetic 
field delayed the onset of wavy instability. So magnetic 
field typically inhibits the oscillatory instability but there 
are situations when magnetic field can stimulate instabil- 
ities [19]. Therefore, a better understanding is necessary 
for this kind of flow. 

Numerical simulations complement experiments be- 
cause very low-Prandtl-number fluids can not be achieved 
in the laboratory. Even the fluids like mercury (P r = 
0.025) do not allow very good visualization in experi- 
ments. However, the simulations are costly in terms of 
computer time. The simulations, due to large number 



of modes, often makes the understanding of the basic 
physics of convective flow difficult. Here, low-dimensional 
models play a very important roll. Low-dimensional 
models are useful for modeling large scale flows. It takes 
much less computer time and yet provide valuable infor- 
mation of the dynamics of the fluid flow. 

In this paper we study the effect of uniform horizontal 
magnetic field applied along a fixed horizontal direction 
(y-axis) using a low dimensional model. In the absence of 
magnetic field the model shows the onset of convection 
in the form of steady two dimensional rolls (2-D). The 
2-D rolls become three dimensional (3-D) via wavy roll 
instability (20[, when the Rayleigh number R is raised 
slightly above critical value R c . The wavy rolls become 
chaotic, if R is raised about 2% above R c for low val- 
ues of Prandtl number P r . These rolls may be aligned 
either along x-axis or along 2/- axis. Two sets of rolls 
compete with each other with further increase in R. The 
horizontal magnetic field does not shift the threshold of 
the primary instability (stationary convection). How- 
ever, it delays the onset of the secondary instability in 
the form of wavy rolls. If Chandrasekhar's number Q, 
which is proportional to the square of the applied hor- 
izontal magnetic field Bo, is raised above certain value, 
the rolls are oriented along the direction of the magnetic 
field in the horizontal plane. The dimensionless num- 
ber R (Q, P r )/(R c QP r ) scales with Q as Q _1 at the on- 
set of wavy instability. If the intensity of the magnetic 
field is sufficient, it suppresses chaotic flow and makes 
the convection periodic for relatively higher values of R. 
A stronger magnetic field suppresses the oscillatory con- 
vection even at higher values of R. The flow becomes 
steady and two dimensional in the presence of stronger 
horizontal magnetic field. The results of the model are 
also compared with those observed in experiments 
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In section II, we describe the physical system together 
with boundary conditions. Section III deals with the 
derivation of the low-dimensional model. The results of 
the model and their comparison with experimental and 
numerical simulations are described in section IV. Con- 
clusions are given in section V. 



II. HYDROMAGNETIC SYSTEM 

As a physical system, a thin horizontal layer of electri- 
cally conducting fluid of thickness <i, uniform kinematic 
viscosity thermal diffusivity ft, magnetic diffusivity A 
and coefficient of volume expansion a is kept between two 
horizontal plates, and is heated uniformly from below. A 
uniform horizontal magnetic field Bq = (0,I?o,0) is ap- 
plied along a fixed direction (say along ?/-axis). Note that 
we have taken x and ?/-axes along horizontal directions 
and z- axis along vertical upward direction. We consider 
the flow of liquid metals which have magnetic Prandtl 
number, P rn (= j) of the order of 10" 5 fl- Therefore 
we set P m = for our study. Now we choose units d 
for length, d 2 jv for time, ^j^- for the induced magnetic 
field and vj3d/ n for temperature, where f3 is the uniform 
temperature gradient between the plates and get the fol- 
lowing set of governing dimensionless equations under 
Boussinesq approximation: 

d t (V 2 v 3 ) = V^ + W^-Qes-fVxVx^b)] 

-e 3 -Vx[(w-V)v-(vV)a;], (1) 
d t uJ3 = V 2 cj 3 + [(lJ'V)vs - (v-V)^] 

+Qe 3 -[Vx(^b)], (2) 

P r [d t 6 + (v • V)0] = v 3 + V 2 #, (3) 

V 2 b = -d y v, (4) 

V-v = Vb = 0, (5) 

where v (x,y,z,t) = (^1,^2,^3) is the velocity field, b 
(x, y, z, t) = (61, 62, bs) the induced magnetic field due to 
convection, (x, z, t) the deviation in temperature field 
from the steady conduction profile, and w = (^1,^2,^3) 
= V x v is the vorticity field in the fluid. The dimension- 
less parameters are: Rayleigh number R = af3gd 4 /vK, 
thermal Prandtl number P r = vj k and Chandrasekhar's 

B 2 d 2 

number Q = — V~ , where q is the acceleration due to 
gravity and po is the reference density of the fluid. £3 
is a unit vector directed vertically upward and V 2 H = 
d X x + d yy is the horizontal Laplacian. We assume the 
horizontal boundaries to be stress-free, thermally and 
electrically good conductors. The boundary conditions 
are then given by 







dvi _ dv2 _ v _ q _ dbi _ db 2 _ 

dz dz ' dz dz 

at z = 0, 1. In what follows we also use one parameter 
r called reduced Rayleigh number, being the ratio of R 
and critical Rayleigh number R c = 277r 4 /4, the threshold 
Rayleigh number for onset of stationary convection. 



III. LOW DIMENSIONAL MODEL 

We apply standard Galerkin technique to derive a sim- 
ple model to describe the convection in presence of hor- 
izontal magnetic field. The spatial dependence of the 
independent fields are expanded in Fourier series com- 
patible with the boundary conditions. The expansions 
are truncated to describe the superposition of two mu- 
tually perpendicular sets of wavy rolls [21] for zero (P r ) 
convection [22[ and extend the model to consider low P r 
convection. The expansions for vertical velocity, vertical 
vorticity and temperature field are as follows: 

^3 = [^ioi(^) cos(kx) + Won(t) cos(ky)] sin(7rz) 
+Wii2(£) cos(kx) cos(ky) sin(27rz) 
+ Win (t) sm(kx) sm(ky) sin(7rz) (6) 

uj 3 = +Z 100 (t) cos(kx) + Z 010 cos(ky) 

+Z\n (t) cos(kx) cos(%) cos(ttz) (7) 
= [T101 (t) cos(kx) + Ton (t) cos(ky)} sin(Trz) 

+ [T n2 (t) cos(kx) cos(%) + T 002 (t)] sin(27rz) 
+^iii (t) sm(kx) sm(ky) sin(7rz) (8) 

The horizontal components of velocity are obtained using 
the solenoidal property of the velocity field. In P m — >• 
limit, the magnetic field fluctuation b is slaved to v, and 
the components of b are determined using Eq. [U We 
then derive the model by projecting the hydrodynami- 
cal equations (PO on the above mentioned modes. The 
model, given in the Appendix, consists of twelve dimen- 
sional coupled first-order nonlinear ordinary differential 
equations for the above Fourier amplitudes. 

IV. RESULTS AND DISCUSSIONS 

We integrate the model using the ode45 solver of MAT- 
LAB to investigate the effect of magnetic field. We set 
k — k c = ir/y/2. The model is integrated in the absence 
of magnetic field using random initial conditions for dif- 
ferent values of the reduced Rayleigh number r = R/R c . 
Once the system attains steady state, magnetic field is 
switched on by setting Q nonzero. We start our discus- 
sion by presenting the results of the model in the absence 
of magnetic field in the next subsection. 



A. Convection in the absence of magnetic field 

In the absence of magnetic field, the model shows the 
onset of convection in the form of steady two dimensional 
rolls as the value of r becomes slightly above unity. As 
r is increased above r G (Q,P r ) = R(°\Q,P r )/ R c , where 
R(°\Q,P r ) is the threshold for oscillatory (secondary) 
instability in the presence of horizontal magnetic field, 
the convective flow becomes three dimensional via oscil- 
latory instability [20|. In the model, the secondary insta- 
bility appears as wavy rolls [23| as r is raised above r Q . 
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With little further increase in r, we find chaotically os- 
cillating wavy rolls either oriented along ?/-axis or along 
x-axis depending upon the choice of initial conditions. 
These wavy rolls are observed for 1.0011 < r < 1.19 for 
P r = .025. For r > 1.19, the model shows chaotic com- 
petition of rolls and squares [2ll-[23| . In the next subsec- 
tions we present the results of the model on the effects of 
magnetic field on the onset of wavy roll instability and 
chaotic flows. 



B. Inhibition of oscillatory instability 

The convection generated magnetic field adds a term 

— Q A -1 (fp-^ only in the dimensionless Navier- Stokes 

equation, where A -1 denotes the inverse of Laplacian. 
This term represents some kind of anisotropic viscos- 
ity [HI . Therefore it inhibits velocity variations along the 
magnetic field direction. As a result, the magnetic field 
inhibits the onset of oscillatory instability and greatly 
enhances the stationary 2-D convection regime. It is ob- 
served that the onset of oscillatory instability depends 
both on P r and Q. Fig. [1(a) clearly shows how the sta- 
bility boundary of the 2-D convection regime pushed to- 
wards higher values of r due to the application of uni- 
form magnetic field in horizontal plan e. This result 
corroborates the experimental [HI, Il6| as well as nu- 
merical [H, G3 results. The threshold for wavy insta- 
bility depends on Prandtl number P r as well as Chan- 
drasekhar's number Q. The combination of dimension- 
less numbers r (Q,P r )/ (QP r ) shows power law behaviour 
(see FigJT(b)), and scales with Q as Q _1 . We have also 



displayed the behaviour of the relative distance from the 
threshold of oscillatory instability (Fiq —R^)/P r in fig- 
ure [2](a) and the angular frequency u Q at the onset of 
oscillatory instability in figure [2(b) as a function of Q. 
Here Rq ^ and are the threshold values of oscillatory 
(secondary) convection in the presence of the external 
magnetic field (Q ^ 0) and in the absence of external 
magnetic field (Q = 0) respectively. The results of our 
model is in good agreement with the numerical [5[ and ex- 
perimental [16[ results showing the power law. The quan- 
tity (R^ - R { o) )/Pr scales with Q as Q a with a = 1.1, 
while the experiments and direct numerical simulations 
gave a = 1.2. The frequency of oscillation at the onset 
of oscillatory instability is linearly proportional to the 
Chandrasekhar's number Q. 




FIG. 2: (a)Scaling of (R { ° } - R { o) )/P r with Q for different 
values of P r . The solid line is parallel to the best fit, and its 
slope is 1.1 ±0.01. (b) Dimensionless frequency oj at the onset 
of oscillatory instability as a function of Q for P r — 0.025. The 
best fit (solid line) shows the linear increase of uo with Q. 
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FIG. 1: (a) Variation of the threshold r (Q, P r ) for wavy insta- 
bility as a function of Prandtl number P r for different values 
of Chandrasekhar's number Q. (b) Scaling of the combination 
r (Q,Pr)/(QPr) with Q for P r = 0.025 (stars) and P r = 0.1 
(circles). The quantity r (Q, P r )/{QP r ) scales as The 
solid line shows the same power law behaviour. 



C. Orientation of convective rolls 

Figure [3] shows the effect of horizontal magnetic field 
on the orientation of the convective rolls. Fourier modes 
Wioi (blue curve) and VFon (green curve) oscillate chaot- 
ically in the absence of the external magnetic field. Fig- 
ure[3](a) shows chaotic wavy rolls along x-axis for r = 1.05 
in the absence of external magnetic field (Q = 0). The set 
of wavy rolls initially along x-axis orient themselves along 
the direction of the magnetic field (see, figure [3fb)) for 
Q = 5. For slightly higher values of the reduced Rayleigh 
number (r = 1.2) , the model shows competition between 
mutually perpendicular sets of wavy rolls (Fig. 13(c) ) in 
the absence of external magnetic field (Q = 0) in stead 
of one set of wavy rolls (Fig. [3(a)). However, the set 
of wavy rolls along x-axis disappears and only the rolls 
along 2/- axis survive (see Fig. [3(d)) for Q = 6. We con- 
clude that the external magnetic field orients the wavy 
rolls along its own direction. 
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FIG. 3: Time series of Wioi (blue curve) and Won (green 
curve) for P = 0.025 showing (a) chaotic wavy rolls along 
x-axis (r = 1.05, Q — 0) and (b) chaotic wavy rolls along 
y-axis (r = 1.05, Q = 5), (c) chaotic competition of mutually 
perpendicular wavy rolls (r = 1.2, Q = 0), and (d) chaotic 
wavy rolls along ?/-axis (r = 1.2, Q = 6). 



FIG. 4: Time series of Wioi and Won at r = 1.05 and P r = 
0.025 for four different values of Q: (a) Chaotic wavy rolls 
along 2/-axis (for Q = 0) are represented by the variation 
of the Fourier mode Wioi (blue). The mode Won (green) 
remains zero as the rolls are already oriented along y— axis, 
(b) periodic wavy rolls (Q = 24), (c) periodic wavy rolls with 
smaller amplitude of oscillation (Q — 38), and (d) stationary 
rolls along y-axis (Q = 45). 



D. Suppression of chaos 



As mentioned above, in the absence of magnetic field 
the model shows three types of chaotic patterns: chaotic 
wavy rolls along ?/-axis, chaotic wavy rolls along x-axis 
and a competition of mutually perpendicular sets of 
chaotic wavy rolls depending upon the choice of initial 
conditions and the value of r. The external uniform hori- 
zontal magnetic field of moderate strength applied along 
2/-axis orients the rolls along its own direction, if it is 
not already oriented along the magnetic field (see, fig- 
ure. [3]). Once the roll-patterns are oriented along the 
field direction, the effect of further increase of magnetic 
field is similar to its effect on a roll-patterns already ori- 
ented along axis. Figure H^a) displays the time series 
of the Fourier modes Wioi (blue) and Won (green) in 
the absence of magnetic field (Q = 0). This represents 
chaotic wavy rolls along 7/-axis. If the external magnetic 
field along 2/- axis is gradually increased to a large value, 
chaotic flow is suppressed and a periodic flow develops. 
The time series of the Fourier mode W101 corresponding 
to the periodic flow is displayed in figure H^b). Further 
increase in the field strength makes the amplitude of os- 
cillation of the mode W101 smaller (see, figure S{c)). The 
flow becomes stationary and two-dimensional (2-D) as 
field strength is increased more. Figure HJ^d) shows the 
time series of W101 corresponding to 2-D stationary so- 
lution. The mechanism of the suppression of chaos and 
reorientation of wavy rolls in the direction of the exter- 
nal magnetic field (y-axis) is displayed very clearly in 
figure [5l It displays the isotherms at z = 0.5 computed 
from the model for different values r and Q. The top 
row shows the various convective patterns in the absence 
of magnetic field (Q = 0): (a) chaotic square patterns 
(r = 1.2), (b) chaotic wavy rolls along y-axis (r = 1.05) 
and (c) chaotic wavy rolls along x-axis (r = 1.05). The 



effect of increasing magnetic field on different patterns 
are shown in three columns of figure [5) As Q increases, 
the wavy rolls orient themselves along the direction of 
the horizontal magnetic field. The chaos is suppressed 
for relatively higher values of Q and the flow becomes 
steady and two dimensional for further higher values of 
Q . The route to chaotic flow from stationary convec- 
tive motion is displayed in figure [6] As Chandrasekhar's 
number Q is decreased from higher to lower value, the 
system shows a series of bifurcations. For much higher 
values of Q the convection is stationary in the form of 
2-D rolls along the direction of the magnetic field. With 
further decrease in Q, there is a Hopf bifurcation from 
stationary state. The convective flow shows time periodic 
behaviour (figure [6ji-b) . As Q is lowered below certain 
value of Q, one more independent frequency develops in 
the flow. Consequently, the convection becomes quasi- 
periodic in time (figure E3>d). Further decrease of Q leads 
to chaotic motion (figure Efe-f). The route to chaos is via 
temporal quasi-periodic convection. The chaotic motion 
becomes more complicated with further decrease in Q 
(see figure [6^-h). 

The amplitude and frequency of oscillatory convection 
due to Hopf bifurcation, as Q is lowered from higher to 
lower values for a fixed value of r, are displayed in fig- 
ure [71 The dimensionless frequency decreases linearly 
with the decrease in Q (figure Et). The square of the am- 
plitude of the oscillatory motion increases linearly with 
the decrease in Q near the onset of Hopf bifurcation (fig- 
ure [TJd). However, the square of the amplitude starts in- 
creasing slower than linear with decrease in Q away from 
the onset of oscillatory behaviour. Ultimately, its growth 
is stopped and it starts decreasing with decreasing Q. 
This may be the influence of quasi-periodic bifurcation 
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FIG. 5: Isotherms at z — 0.5 for different values of r and 
Q. (a) The left column displays isotherms for r = 1.2, and 
for three values of Q. The chaotic competition of rolls and 
squares (top row) for Q = 0, time-periodic wavy rolls (mid- 
dle row) along y-axis for Q = 60, and stationary 2-D rolls 
(bottom row) for Q = 80. (b) The middle column displays 
isotherms for r = 1.05 and three different values of Q. Chaotic 
wavy rolls (top) along 2/- axis for Q — become time-periodic 
wavy rolls (middle) for Q — 24 and finally spatially periodic 
stationary 2-D rolls (bottom) along y-axis for Q = 45. (c) 
Isotherms for r = 1.05 but with initial conditions different 
than used in (b) are displayed in right column. Chaotic wavy 
rolls along x-axis (top) orient themselves along y-axis as time- 
periodic wavy rolls (middle) for Q — 38, and finally become 
2-D stationary rolls along y-axis for Q = 45. 



present in the model. 

We have also plotted the Nusselt number obtained 
from the model as a function of Rayleigh number R in 
figure [51 The time averaged Nusselt number at the onset 
of convection has strong similarity with the behaviour of 
Nusselt number computed in the direct numerical sim- 
ulations [6]. In addition, the results of the model also 
show the convergence of Nusselt numbers for higher val- 
ues of R as observed in simulation [6] for various values 
of Q. For intermediate values of R, the Nusselt number 
obtained from the model shows differences with the same 
obtained in DNS. This may be due to severe truncation 
of the Fourier expansion of various fields in the model. 
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FIG. 6: Quasi-periodic route to chaos as a function of Q for 
P r = 0.025 and r = 1.05. The left column displays projection 
of the phase space trajectories on Wioi — ^oio plane for four 
different values of Q. The right column shows power spectral 
density (PSD) of the Fourier mode Wioi as a function of the 
dimensionless frequency / for the corresponding values of Q. 
Time-periodic behaviour of the convective flow (a) for Q = 24 
and one peak in PSD corresponding to the periodic flow (b). 
The temporal quasi-periodic behaviour of the flow (c) for Q = 
22.8 and two peaks (d) confirming the quasi-periodic flow. 
Chaotic flow (e) for Q = 22.4 and its PSD (f). The chaotic 
flow becomes more complex for Q = 20 (g) with flattening of 
peaks in PSD (h). 




FIG. 7: (a) Nondimensional frequency / of the wavy convec- 
tion with Q at a fixed value of the reduced Rayleigh number 
r = 1.05. The best fit (solid line) shows the linear dependence 
of / with Q. (b) Square of the amplitude of the periodic wavy 
convection as a function of Q at r = 1.05. 



V. CONCLUSIONS 

In this paper we presented a low-dimensional model 
to study the effect of uniform horizontal magnetic field 
on Rayleigh-Benard convection in low- Prandtl- number 
fluids. The model is capable of capturing convective 
patterns in the form of two sets of mutually perpendic- 
ular straight (2-D) rolls as well wavy (3-D) rolls. The 
external field affects the flow patterns significantly. It 



inhibits the onset of wavy rolls, which appear at the 
secondary instability. The stationary straight rolls at 
primary instability as well as wavy rolls at secondary 
instability orient themselves along the direction of the 
magnetic field, if it is sufficiently strong. The threshold 
value of the Rayleigh number at the onset of wavy in- 
stability shows power law behaviour. Stronger magnetic 
field causes the suppression of chaotic flow. The results 
obtained from this simple model have good qualitative 
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agreement with those observed in experiments and 
simulations. 



FIG. 8: Plot of convective heat flux (Nu — 1) as a function 

of relative Rayleigh number (R - R c ) for P r = 0.025 (a) and pp thanks M . K . Verma, IIT Kanpur, India, for useful 

P r = 0.1 (b) for several values of Q. discussions. 
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Appendix : The model 



u = 


--Wl- 




97T 4 7T 

f- 2QA)U + — rT + -V 2 MXJ 




+^-(ttV 1 C-2X)Z, 


Vi = 


-> 2+ 


27 

Q)Vi + —7r 4 rSi - CU • Z, 

O 


v 2 = 




+ Q)^ 2 + ^7r 4 r,S 2 -y'7 1 C/ 2 , 


z = 


-\(«H + 


2 

o o 


X = 




Q)X + ttU-Z, 


f = 


-2^T + 

2P r 


yU+^SiZ+ J(5 2 B + 4IF)U, 


Si = 




- ^Vi - CZ • T + nViY, 


S2 = 




.i.y 2 _|u-BT, 


Y = 


Pr 


7T 7T 

-U-T--^!, 



where U = {U U U 2 ) T = {W WI , W mi ) T , (V U V 2 ) = 
(Wiu,Wu 2 ), Z = (Z U Z 2 ) T = (Z 010 ,Z W0 ) T , T = 
(Ti,T 2 ) T = (Tioi,Ton) T , (Si,S 2 ) = (Tin, Tm), X = 

Zm, Y = T 002 , A=(J5)' B= (io)' C = 

(o 1 -°i)' D= (J u)' andl = (J !)■ Thedot 

(•) operation implies standard inner product, and the su- 
perscript T denotes Transpose of a matrix. 



